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Abstract 

Background: Tlie main cause of deatli of breast cancer patients is not tine primary tumor itself but the metastatic 
disease. Identif/ing breast cancer-specific signatures for metastasis and learning more about the nature of the genes 
involved in the metastatic process would 1) improve our understanding of the mechanisms of cancer progression 
and 2) reveal new therapeutic targets. Previous studies showed that the transcriptional regulator megakaryoblastic 
leukemia-1 (Mkll) induces tenascin-C expression in normal and transformed mammary epithelial cells. Tenascin-C is 
known to be expressed in metastatic niches, is highly induced in cancer stroma and promotes breast cancer 
metastasis to the lung. 

Methods: Using HCl 1 mammary epithelial cells overexpressing different Mkll constructs, we devised a subtractive 
transcript profiling screen to identify the mechanism by which Mkll induces a gene set co-regulated with 
tenascin-C. We performed computational analysis of the Mkll target genes and used cell biological experiments to 
confirm the effect of these gene products on cell behavior. To analyze whether this gene set is prognostic of 
accelerated cancer progression in human patients, we used the bioinformatics tool GOBO that allowed us to 
investigate a large breast tumor data set linked to patient data. 

Results: We discovered a breast cancer-specific set of genes including tenascin-C, which is regulated by Mkll in a 
SAP domain-dependent, serum response factor-independent manner and is strongly implicated in cell proliferation, 
cell motility and cancer. Downregulation of this set of transcripts by overexpression of Mkll lacking the SAP domain 
inhibited cell growth and cell migration. Many of these genes are direct Mkll targets since their promoter-reporter 
constructs were induced by Mkll in a SAP domain-dependent manner. Transcripts, most strongly reduced in the 
absence of the SAP domain were mechanoresponsive. Finally, expression of this gene set is associated with 
high-proliferative poor-outcome classes in human breast cancer and a strongly reduced survival rate for patients 
independent of tumor grade. 

Conclusions: This study highlights a crucial role for the transcriptional regulator Mkll and its SAP domain during 
breast cancer progression. We identified a novel gene set that correlates with bad prognosis and thus may help in 
deciding the rigor of therapy. 
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Background 

Most breast cancer patients die from tumor metastases 
and not from the primary tumor itself. Thus, the identi- 
fication of genes and signaUng pathways influencing the 
metastatic process are of utmost importance. Once the 
mechanisms leading to metastasis are uncovered, they 
can in the future serve as a rational basis for prognosis 
and intervention. From the beginning of its discovery, 
tenascin-C has been strongly associated with tumorigen- 
esis and cancer progression in many different types of 
tumors (reviewed in [1,2]). Tenascin-C was not only 
enriched in breast cancer tissue [3,4], but its high ex- 
pression was part of a gene signature of breast cancers 
metastasizing to the lung [5]. There is strong evidence 
that tenascin-C contributes to the metastatic behavior of 
breast cancer cells [6] by providing a niche for their 
settlement in the lung [7,8]. The source of tenascin-C 
can be the tumor cells themselves as well as the stromal 
cells of the cancer microenvironment. Downregulation 
of tenascin-C by miR-335 or shRNA in human cancer 
cells in a mouse xenograft model inhibits metastasis for- 
mation [7], and in tenascin-C-deficient mice, metastasis 
formation of tenascin-C positive cancer cells is also 
suppressed [9]. 

There are many signaling pathways inducing tenascin-C 
expression (reviewed in [10]). Among these, mechanical 
strain appUcation in vivo as well as to cells in culture is 
a potent stimulus to induce tenascin-C expression in fi- 
broblasts [11,12]. We have recently shown that induc- 
tion of tenascin-C by cyclic mechanical strain requires 
the action of Midi [13]. Midi is a member of the 
myocardin-related transcription factor family (MRTF) 
and a well-known transcriptional co-activator of serum 
response factor (SRF) [14-16]. SRF target genes, which 
are regulated upon recruitment of MRTF cofactors, en- 
code proteins involved in actin cytoskeletal function 
that can either be structural (for example, actin) or re- 
lated to actin dynamics (for example, talin 1) (reviewed 
in [17,18]). However, Midi -mediated stretch-induced 
tenascin-C expression in fibroblasts did not require 
SRF, but instead depended on the potential DNA- 
binding SAP domain of Midi. This implies a novel 
mode of Midi action as a bona fide transcription factor 
in mechanotransduction [13]. Interestingly, normal and 
transformed mouse mammary epithelial cells also ap- 
pear to be highly sensitive to Midi signaling, respond- 
ing to Midi overexpression with several fold induction 
of tenascin-C [13]. 

The present study was designed to find SAP-dependent 
Midi target genes co-regulated with tenascin-C and to 
analyze whether such genes could be indicative of specific 
physiological states of cells that might be controlled by 
mechanotransduction. For our study, we made use of the 
HCll mammary epithelial cell line. HCll cells are capable 



of both self-renewal and differentiation and can be cul- 
tured for unlimited time in an undifferentiated state [19], 
the condition we used in our study. HCll cells can recon- 
stitute the ductal epithelium of a cleared mammary fat 
pad in vivo with ductal, alveolar and myoepithelial cells, 
illustrating their stem cell abilities [19,20]. In addition, 
HCll cells contain a mutated p53 gene that not only in- 
creases the replicative potential of stem cells but confers 
predisposition to mammary carcinoma [21]. Undifferen- 
tiated HCll cells share transcriptome signatures with 
human breast cancer [22], supporting the relevance of 
this model for breast cancer-related studies. We there- 
fore concluded our study by investigating whether the 
genes co-regulated with tenascin-C would also be impli- 
cated in breast cancer progression. 

Results 

Screen for SAP-dependent Mkll target genes 

We devised a screening method to identify genes co- 
regulated with tenascin-C by Midi in a SAP domain- 
dependent manner without involvement of SRF. For this 
purpose, we used HCll mammary epithelial cells that 
react strongly to the overexpression of Midi with in- 
duction of tenascin-C expression [13]. We compared 
three HCll strains that either overexpress the C-terminal 
red fluorescent protein (RFP) -tagged full length Mkll 
(HCll-FL), Mldl-RFP with a mutated SRF-interaction site 
(HCll-mutBl) or Mkll-RFP with a deletion of the SAP 
domain (HCll -ASAP). None of the three Midi variants 
appear to be toxic to the cells, as we have not observed 
any changes in viability or cell morphology. HCll-FL cells 
were shown to overexpress Mkll 7.1-fold above the en- 
dogenous Midi present in parental HCll cells [13], and 
were used as control cells in our study. All cell strains 
were FACS sorted to express similar levels of Mldl-RFP 
proteins. These cells were used for transcript profiling 
and gene lists of interest were established as shown in 
Figure lA, B. A scatter plot (Figure lA) of all transcripts 
expressed in HCll-mutBl versus HCll-FL control cells 
(y-axis) and all transcripts expressed in HCll- ASAP 
versus HCll-FL control cells (x-axis) shows that a large 
majority of transcripts does not differ significantly be- 
tween the three cell strains (log fold change (FC) ~ 0; 
black dots). Setting the threshold to a 2-fold reduction 
(logFC = -1; grey lines), three gene sets can be distin- 
guished: 1) blue dots represent genes that are lower in 
HCll-mutBl than in HCll-FL control cells, but are 
unaffected in HCll-ASAP cells, thus representing typ- 
ical SRF/Mldl target genes; 2) green dots represent 
genes that are lower in HCll-ASAP than in HCll-FL 
control cells, but are unaffected in HCll-mutBl cells 
(this gene set includes tenascin-C); and 3) red dots indi- 
cate genes with reduced expression in both HCll-mutBl 
and HCll-ASAP cells compared to HCll-FL control cells. 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 Screen for SAP-dependent Mkll target genes and their implication in cancer. (A) Scatter plot and (B) Venn diagram representing 
classification of Mkll target genes into three groups: SRF-dependent/SAP-independent (blue), SRF-dependent/SAP-dependent (red) and 
SRF-independent/SAP-dependent (green). The scatter plot (A) represents the log fold change (logFC) in gene expression in HCl 1-ASAP versus 
HCl 1-FL control cells (x-axis; ASAP vs. FL) and between HCl l-mutBl versus HCl 1-FL control cells (y-axis; mutBl vs. FL). Each dot represents a 
probeset, and the one for tenascin-C is highlighted (Tnc). The vertical and horizontal lines in the chart denote the 2-fold change cutoff (logFC 
= -1). The Venn diagram (B) represents the number of probesets for transcripts, which are more than 2-fold reduced in either HCl 1 -mutBl or 
HCl 1-ASAP cells when compared to HCl 1-FL control cells. Boxes below the Venn diagram indicate the cell strains that have reduced levels of 
the respective transcripts. (C D) Functional analysis for the three Mkll -regulated gene sets performed using the IPA software. The high-level 
functional (C) and disease (D) categories are displayed along the x-axis of each bar chart. The y-axis displays the -log of the P-value determined 
by right-tailed Fisher's exact test. The P-value is a measure of the likelihood that the association between a set of genes in each dataset and a 
related function or disease is due to random association. The grey vertical line denotes the cutoff for significance (P = 0.05; -logP = 1 .3). 



Thus, this approach enabled us to form three gene sets 
that were distinct from the large majority of transcripts 
and were dependent for expression on the Bl site of 
Mkll, the SAP domain, or both. The three groups pre- 
sented by a Venn diagram (Figure IB) contain 141 pro- 
besets for transcripts that depended on the function of 
the Bl site but not the SAP domain for their induction, 
113 probesets for transcripts that depended on both of 
these Mkll domains and a third group of 205 probesets 
for transcripts co-regulated with tenascin-C that did not 
require an interaction of Mkll with SRF but depended 
on the SAP domain for induction (complete probeset 
Usts and annotations are found in Additional file 1: 
Table SI, Additional file 2: Table S2 and Additional file 3: 
Table S3). This analysis revealed that the SAP-dependent 
mechanism of tenascin-C regulation by Mkll is shared 
by a large cohort of genes. Below the Venn diagram, we 
indicated which cells were deficient in the respective 
transcripts. Thus, the typical SRF/Mkll target genes are re- 
duced in HCll-mutBl cells, while the SRF-independent/ 
SAP-dependent genes are reduced in HCl 1-ASAP cells. 
The intermediate group that requires both Mkll activities 
is reduced in both the HCll-mutBl and HCll-ASAP 
cells. 



The SAP-dependent Mkl1 target genes are implicated in 
cancer 

Functional analysis of the three gene lists using the 
IPA software revealed different molecular and cellular 
functions (Figure IC) and different disease associations 
(Figure ID) for the three types of gene signatures. Thus, 
the SRF-dependent/SAP-independent signature implicated 
a function of these genes in cellular movement and the 
linked diseases included connective tissue disorders, in- 
flammatory disease and skeletal and muscle disorders, 
which are the main features known to be regulated by 
SRF/Mkll interaction [23-25]. The SRF-dependent/SAP- 
dependent group of genes includes as major functions 
post-translational modification, protein degradation and 
protein synthesis, and the top disease association is cancer. 
Finally, the genes of the SRF-independent/SAP-dependent 



group were associated with extremely high significance 
with cell cycle and cancer (-logP > 25 and > 30, respect- 
ively), while the SRF/Mkll target genes were associated 
with the same two categories at low significance only 
(-logP > 2 and > 7, respectively). These data imply that 
SAP-dependent induction of transcription by Mkll may 
counteract the known differentiation-promoting effect of 
SRF/Mkll -induced transcription. A list of SAP-dependent 
genes with published cancer-related functions, whose 
transcripts were downregulated more than 3-fold in 
HCll-ASAP compared to HCl 1-FL control cells, is pre- 
sented in Table 1. To confirm that these transcripts are 
indeed differentially expressed in the different HCll 
cell strains, qRT-PCR analysis was performed using 
cDNA from three different batches of the respective 
HCll strains. Differences in gene expression between 
HCll-ASAP and control cells are presented in Table 1 
and in more detail in Additional file 4: Figure SI. The 
qRT-PCR results agreed with the data obtained by tran- 
script profiling. We also tested the SAP-dependent gene 
expression in the HCll strains when grown in the pres- 
ence of serum. It is interesting to note that in the presence 
of 3% FCS, these transcripts remained strongly reduced in 
HCll-ASAP compared to control cells (Table 1). Thus, 
the induction of these genes seems to depend mainly on 
whether the SAP domain is present in the transfected 
Mkll construct. 

In addition, we monitored changes in the expression 
of some of the SRF-independent/SAP-dependent Mkll 
targets on a protein level. In agreement with the changes 
seen at the transcript level, we confirmed the reduction 
of tenascin-C, Wispl and Nox4 proteins in cells overex- 
pressing the ASAP-Mkll construct compared to the 
HCll-FL control and HCll-mutBl cells (Additional file 4: 
Figure S2). Using zymography, we found that Mmp2, a 
gene that was not affected by Mkll overexpression at 
the transcript level was highly expressed in all three cell 
strains, whereas Mmp3 and/or 12, which belonged to 
the SRF-dependent/SAP-dependent gene set, were al- 
most completely lacking in HCll-mutBl as well as 
HCll-ASAP cells, corresponding to the data obtained 
by transcript profiling. 
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Table 1 SAP-dependent Mkll target genes 



Gene Description Fold Reduction in Functions 

HC11-ASAP vs. HC11-FL cells 
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in 0.03% FCS 


qRT-PCR in 
0.03% FCS 


qRT-PCR 
in 3% FCS 
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Cell cycle regulation [26], cell motility and cancer 
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Cell growth, differentiation and migration [29], tumor 
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3.63*** 


5.70*** 


14.84** 


Cell growth and motility [31], role in arthritis [32] and cancer 
[31] 


r\l ID 


IxcldLIII J, II 1 Lcl 1 1 IcUldLc illdnicML 

protein 


3.73*** 


2.74*** 


8.02*** 


Protein synthesis, epithelial cell growth and 
differentiation [33,34] 




7R1 041 7H1 ^Rik 
zo 1 1 / n 1 or\ii\, 

PCNA-associated factor 


3.91*** 


1 .89*** 




DNA repair and cell cycle regulation, cell survival and 
proliferation, tumorigenesis [35-37], hematopoiesis [38] 


Assi 


Argininosuccinate synthetase 1 


4.23*** 


3.89** 


2.72** 


Regulation of nitric oxide production and cell viability [39,40] 


Cd34 


CD34 antigen, stem cell antigen 


4.ZJ 


1 U.D 1 


1 ./Z 


Vessel development and function [41], tumor growth [42,43] 


Wispl 


WNTl inducible signaling 
pathway protein 1, ECM protein 


1 


Z.J4 


4.UD 


Cell proliferation and survival, ECM deposition and turnover, 
EMT, tumorigenesis, tissue remodeling [44] 


Mcm6 


Minichromosome maintenance 
complex component 6 


4.4z 


Z.OD 


1 .oU 


Cell cycle regulation [45] 


Carl 2 


Carbonic anyhydrase 12 
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Cell survival under hypoxic conditions, tumor-associated cell 
migration and invasion [46,47] 
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Hyaluronectin, TIP30, 
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245 27*** 
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metastasis [49] 


Kif26b 


Kinesin family member 26B 


6.33*** 


8.36*** 


61.22*** 


Regulation of adhesion and cell polarity in kidney 
development [50] 


SRF-dependent genes 










Lox 


Lysyl oxidase, ECM protein 


4.61*** 


4.70** 


12.04*** 


ECM turnover, connective tissue remodeling and repair, 
tumor progression and metastasis [51,52] 


Mmpl2 


Matrix metallopeptidase 12, 
metalloelastase 


12.01*** 


23.49*** 


4.90** 


ECM degradation in tissue remodeling [49] and 
tumorigenesis [53] 


MmpS 


Matrix metallopeptidase 3, 
stromelysin-1 


15.64*** 


14.70*** 


2.08** 


ECM degradation in tissue remodeling [49] and 
tumorigenesis [54] 



Abbreviations: ECM extracellular matrix protein, PCNA proliferating cell nuclear antigen, EMT epithelial-to-mesenchymal transition. 
***P < 0.001, **P < 0.01, Student's f test. 



SRF-lndependent/SAP-dependent transcripts represent 
direct Mkll target genes 

Since we have previously shown that the SAP domain of 
Midi interacts with the proximal promoter of tenascin-C 
to induce its transcription [13], we tested whether this was 
also the case for other transcripts of the same group. 
The promoters of the SRF-independent/SAP-dependent 
genes listed in Table 1 encompassing at least 500 bp up- 
stream of the transcription start site (TSS) were fused 
to the secreted alkaline phosphatase (SEAP) reporter 
gene of pSEAP2-Basic. We tested the induction of each 
promoter-reporter construct by co-transfection with FL- 
Mkll (Figure 2A). This revealed that the majority of the 
new promoters tested (8 out of 12) were induced at least 
2-fold by Mkll in comparison to co-transfection with 



an inactive Mkll devoid of the transactivation domain, 
indicating that these are indeed direct Mkll target 
genes. The promoter constructs that did not respond to 
Mkll overexpression may represent genes that are in- 
directly regulated by Mkll, or the relevant promoter 
regions were not contained in the constructs tested. 
Next, we investigated whether the induction was SAP- 
dependent and SRF-independent by comparing the re- 
porter activity after co-transfection with mutBl- versus 
ASAP-Mkll variants (Figure 2B). Indeed, the promoter- 
reporter constructs induced by FL-Mkll were also 
strongly induced by mutBl-Mkll, but not by ASAP-Mkll. 
In contrast, the promoter construct for Acta2, a gene from 
the SRF-dependent/SAP-independent gene set was strongly 
induced by ASAP-Mkll but not by mutBl-Mkll, as 
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Figure 2 SRF-independent/SAP-dependent transcripts represent 
direct Midi target genes requiring the SAP domain of Midi to 
induce transcription from their proximal promoter. (A) The 

indicated promoter constructs tliat contained at least 500 bp 
upstream of tine transcription start site (JSS) and were linl<ed to tlie 
secreted all<aline pliospliatase (SEAP) reporter gene, were 
cotransfected in HCl 1 cells together with an inactive Mkll devoid of 
the transactivation domain [13] or the FL-Mkll construct. SEAP activity 
is expressed as fold induction above the level obtained with the 
inactive Mkll. In addition to Tnc, for 8 out of the 12 new promoters 
tested, induction greater than 2-fold (indicated by the red line) was 
obtained. Values are means ± SEM from three to seven independent 
experiments. (B) HCl 1 cells were cotransfected with the indicated 
promoter constructs that were either > 500 bp or shortened to 200 bp 
upstream of theTSS, and with vectors encoding the indicated mutant 
Mkll constructs. SEAP activity is normalized and expressed as in 
(A). Means ± SEM from at least three independent experiments and 
significant differences between either mutBl- and ASAP-Mkll- 
transactivated promoter constructs or between the longer and shorter 
promoter constructs transactivated by mutBl -Mkll, ***P < 0.001, 
**P< 0.01, *P< 0.05 are shown. 



expected for a typical SRF/Mkll target gene [18,55,56]. All 
promoters that revealed SAP-dependency were shortened 
to 200 bp upstream of the TSS to test whether this was 
sufficient to relay the Mkll response, as it has been seen 
previously for tenascin-C [13]. With the exception of Krt5 



and Nox4, for which some activity was lost by shortening 
the promoters, the 200 bp proximal promoters of all other 
genes tested were induced equally well as the longer con- 
structs (Figure 2B). Thus, we conclude that there are many 
genes that are regulated similarly as tenascin-C requiring 
the SAP domain of Mkll to induce transcription from 
their proximal promoter. 

The different HCl 1 cell strains proliferate at different 
rates and show distinct migration behaviors 

Next, we tested whether the differential gene expression 
seen in the different HCll strains overexpressing either 
FL-, mutBl- or ASAP-Mkll constructs have functional 
consequences on their behavior. Since most of the SAP- 
dependent transcripts are proposed to have a function 
in cancer, we decided to analyze two main functions im- 
portant for cancer progression: proliferation and cell mi- 
gration. An approximately equal overexpression of the 
different Mkll protein variants in the HCll cell lines was 
confirmed by Western blot analysis (Figure 3A). An HCll 
cell strain stably transfected with an empty vector [13] 
expressing only endogenous Mkll (below the detection 
limit in Figure 3A) was also included in these studies. 
The proliferation rates of the HCll strains were ana- 
lyzed using a 5-bromo-2'-deoxyuridine (BrdU) incorp- 
oration assay. The incorporated BrdU was measured 
immediately after plating (0 h) as well as at 24, 48, 72 
and 96 h. Compared to empty vector-, FL- or mutBl- 
transfected HCll strains, there was a significant decrease 
in BrdU uptake into newly synthesized DNA in HCll- 
ASAP cells over the entire time period tested (Figure 3B). 
To investigate cell motility, we used a transfilter migration 
assay. Similarly to the effect on cellular proliferation, the 
expression of ASAP-Mkll significantly inhibited HCll cell 
migration by 2.7-fold compared to endogenous or full 
length Mkll expression, and more than 3.5-fold compared 
to mutBl-Mkll expression (Figure 3C). 

Thus, overexpression of FL-Mkll protein in HCll 
cells did not affect their behavior. However, overexpres- 
sion of ASAP-Mkll led to a significant reduction in the 
proliferative and migratory ability of HCll epithelial cells, 
either through a dominant-negative effect of ASAP-Mkll 
on SRF-mediated action and/or a positive impact of the 
SAP-dependent Mkll target genes on these functions 
important for cancer progression. 

SAP-dependent Mkll target genes are 
mechanoresponslve 

We have previously found that the SAP-dependent in- 
duction of tenascin-C was triggered by applying mech- 
anical strain to fibroblasts. Mammary epithelial cells are 
also exposed to mechanical strains, both during normal 
development, pregnancy and lactation, as well as under 
pathological conditions such as in cancer. Therefore, we 
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Figure 3 The different HC1 1 cell strains proliferate at different 
rates and show distinct migration behaviors. (A) Immunoblot 
with mAb65F13 of Mkll proteins in wliole-cell extracts from tine 
empty vector-, FL-, mutBI - or ASAP-transfected HCl 1 strains. 
Anti-Gapdii served as loading control. Endogenous Mkll protein was 
below the detection limit in empty vector cells. (B) SAP-dependent 
proliferation of HCl 1 mammary epithelial cells. Proliferation rates of 
the four HCl 1 cell strains were assessed by BrdU incorporation into 
newly synthesized DNA immediately after plating (0 h) as well as at 
24, 48, 72 and 96 h. Means ± SD from three independent 
experiments and significant differences to the HCl 1-ASAP cells, 
^*^P < 0.001, ^^P < 0.01, *P < 0.05 are shown. (C) SAP-dependent 
migration of HCl 1 mammary epithelial cells. Cell migration of the 
four HCl 1 strains was evaluated by Transwell migration assay using 
filters with 8 \im pore size. Quantification of the cell migration was 
measured by the area on the lower side of the filter covered with 
cells. Above the bar graph, a photo of fixed and stained cells seeded 
in parallel in a 24-well plate is shown as a seeding control, and 
representative photos of fixed and stained cells of each of the cell 
strains that have migrated to the lower side of the filter, are shown 
below (bar, 200 \im). Data and statistical significance are expressed 
as in (B). 



tested whether tenascin-C and other members of the 
SAP-dependent Mkll -induced gene set are mechanore- 
sponsive in HCll cells. We tested two paradigms: 1) static 
strain that was shown to induce c-fos, a very prominent 
mechanoresponsive gene in HCll cells [57] that we used 
as a control, and 2) cyclic strain. While we were able to 
confirm induction of c-fos by applying static strain at 20% 
for 1 h, there was no induction of tenascin-C under these 
conditions compared to cells at rest (Figure 4A). However, 
using 15% cyclic strain at a frequency of 0.3 Hz for 1 h, we 
found that not only the control gene c-fos but 11 out of 
16 SAP-dependent genes, including tenascin-C were sig- 
nificantly upregulated above the expression levels obtained 
in resting cells (Figure 4). Even though significant, the in- 
duction of tenascin-C was minimal (Figure 4A) compared 
to 18-fold upregulation for Adamtsl6 or 10-fold upreg- 
ulation for Lox (Figure 4B), both of which are enzymes 
involved in extracellular matrix (ECM) remodeling and 
cancer progression [31,58]. Being mechanoresponsive, 
the SAP-dependent Mkll target genes might be activated 
in stiff tumor tissue, which further confirms their relation 
with cancer. 

The SRF-independent/SAP-dependent genes represent a 
bad prognostic signature for breast cancer patients 

In order to investigate whether the SRF-independent/SAP- 
dependent genes were prognostic of accelerated cancer pro- 
gression in human patients, we used the bioinformatics tool 
Gene expression-based Outcome for Breast cancer Online 
(GOBO) that allowed us to investigate a breast tumor 
data set containing 1881 samples analyzed by Affyme- 
trix Human Genome U133A arrays. GOBO is designed 
to assess gene expression levels and association with out- 
come of single genes or gene sets in multiple subgroups of 
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Figure 4 SAP-dependent Mkll target genes are 
mechanoresponsive. (A) Effect of static (20%) and cyclic (15%, 
0.3 Hz) strain on Tnc and c-fos mRNA levels. HCl 1 cells were 
cultured on either growth factor reduced matrigel matrix- or 
fibronectin-coated silicone membranes in 0.03% serum-containing 
medium for 24 h before applying static or cyclic strain for 1 h. Cells 
cultured under the same conditions and not exposed to mechanical 
stimulation were used as a resting control. The two types of coating 
gave identical results under the indicated experimental conditions. 
Total RNA was extracted and qRT-PCR was performed for Tnc and 
c-fos mRNA levels. Values normalized to Gapdh are expressed 
relative to the values of resting cells. Data represent means ±SD 
from three independent experiments. Significant differences to the 
resting control, ***P < 0.001, ^^P < 0.01, *P < 0.05. (B) SAP-dependent 
genes respond to cyclic strain. HCl 1 cells were stretched and mRNA 
analyses were performed as described in (A). Data and statistical 
significance are expressed as in (A). 



this breast cancer data set [59]. Here, we analyzed two 
sets of genes, namely the SRF/Mkll -induced gene set 
(SRF-dependent/SAP-independent) and the SAP-dependent 
gene set (SRF-independent/SAP-dependent) containing 
tenascin-C. The analysis was performed across tumor 
samples stratified according to PAM50 subtypes [60], 
estrogen receptor (ER) -status and histological grade. In 
contrast to the SRF/Mkll target genes that were predom- 
inantiy associated with tumors classified as normal-like 



and with lower histological grades (1 and 2) (Figure 5A), 
elevated expression of SAP-dependent genes was associ- 
ated with extremely high significance (P < 0.00001) with 
typical high-proliferative poor outcome classes in breast 
cancer, such as basal-like, HER2-enriched, luminal B, ER- 
negative and histological grade 3 tumors (Figure 5B). Next, 
a functional correlation analysis to find a possible inter- 
connection between the SAP-dependent Mkll target 
genes was performed using the GOBO tool (Additional 
file 4: Figure S3). This analysis explores the correlation 
of expression of individual genes in our gene sets with 
eight different co-expressed gene modules emulating 
breast cancer-specific as well as general tumor biological 
processes [61]. Interestingly, whereas the gene set of SRF/ 
Mkll targets did not show a significant correlation with 
any of these modules, the genes in the SAP-dependent 
gene set were correlated with a very high significance 
(P < 0.00001) with two proliferation modules - mitotic 
checkpoint and mitotic progression. Both modules con- 
tain genes related to central mitotic processes involved 
in either the regulation of the M-phase and the mitotic 
checkpoint or in carrying out the M-phase. Finally, the 
association of our gene sets with outcome using distant 
metastasis free survival (DMFS) as an endpoint and 10- 
year censoring was analyzed. The survival analysis was 
performed in all tumors for which DMFS follow-up is 
available (1379 cases), as well as in 21 groups that were 
stratified based on gene expression subtypes (PAM50 
classifier), ER-status, lymph node (LN) -status, histological 
grade, and treatment status. Samples in the whole cancer 
data set (1881 patients) were stratified into three quan- 
tiles, low, intermediate and high, based on SRF-dependent/ 
SAP-independent or SRF-independent/SAP-dependent gene 
expression. Interestingly, high expression of SRF/Mkll - 
induced genes was associated with a better clinical out- 
come for all tumors, as well as for LN-negative and 
untreated tumors compared to low and intermediate ex- 
pression of these genes (Figure 6A). In contrast, both high 
and intermediate expression of the SAP-dependent genes 
was associated with bad clinical outcome in all tumors, 
and particularly in LN-negative, systemically untreated, 
ER-positive, Grade 1 and 2 tumors (Figure 6B). Similar re- 
sults were obtained for the typical breast cancer gene 
CCNBl by Ringner et al. [59]. The Kaplan-Meier survival 
analyses were supported by the corresponding multivariate 
analyses (Figure 7 A, B). The hazard ratio for the variate 
Grade shows statistical significance, proving that the in- 
fluence of high SAP-dependent gene expression on pa- 
tient survival is independent of tumor grade. Among all 
tumors for which DMFS data are available, a hazard ra- 
tio of 0.44 (95% CI = 0.28-0.68; P = 0.0003) for the low 
SRF-independent/SAP-dependent tercile was detected 
compared to the high SRF-independent/SAP-dependent 
tercile (Figure 7B, all tumors). This indicates that patients 
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Figure 5 SAP-dependent Mkll target genes are associated witii typical high-proliferative poor outcome classes in breast cancer. The 

expression levels for the SRF-dependent/SAP-independent (A) and SRF-independent/SAP-dependent (B) gene sets are analyzed across the 
1881 -sample breast cancer data set stratified according to PAM50 subtypes (left panels), estrogen receptor (ER)-status (middle panels) and 
histological grade (right panels), and represented by box plots using the GOBO bioinformatics tool. The number of tumors in each breast cancer 
subtype and the significant difference in gene expression (P-value calculated using ANOVA) between them are shown above the box plots. 



with tumors expressing high levels of the SAP-dependent 
genes are more than twice as likely to develop metastatic 
disease. Similar hazard ratios, in the range of 0.28-0.44 
for the low tercile compared to the high tercile were also 
detected among subgroups of untreated, LN-negative, 
ER-positive, Grade 1 and 2 tumors (Figure 7B). Thus, 
the association of high SRF-independent/SAP-dependent 
gene expression with reduced DMFS among patients not 
receiving adjuvant therapy, as well as among LN-negative, 
ER-positive, Grade 1 and 2 patients indicates that in- 
creased expression of the SAP-dependent Mkll target 
genes plays a significant role in the natural metastatic 
progression of non-aggressive towards highly aggressive 
breast cancer in human patients. 



Discussion 

Given the heterogeneity of mutations in tumor cells, it 
becomes increasingly clear that not only individual 
genes but pathways govern the course of tumorigenesis 
and cancer progression [62]. We have recently shown 
that induction of tenascin-C by cyclic mechanical strain 
required the action of the potential DNA-binding SAP 
domain of Mkll independently of an interaction of 
Mkll with SRF [13]. Now, we report a screen for genes 
co-regulated with tenascin-C by the same SAP-dependent 
and SRF-independent mechanism in mammary epithelial 
cells. This screen reveals a set of SAP domain-dependent 



Mkll target genes with a strong implication in cell prolif- 
eration, cell motility and cancer. 

To date only a few studies have shown that Mkll is 
implicated in cancer-related processes (reviewed in [63]) 
and most of them have concentrated on the SRF/Mkll 
signaling for the induction of individual genes [64-67]. 
The first study reporting that depletion of Mkll/2 proteins 
reduced motility, invasion and colonization of metastatic 
tumor cells in an experimental in vivo metastasis assay 
[64] was further supported by the discovery of the Mkll- 
binding protein, suppressor of cancer cell invasion (SCAT), 
which inhibited SRF/Mkll-mediated expression of |3l in- 
tegrin [68] . Since then, several studies describing opposing 
biological effects for Mkll appeared. For instance, several 
antiproliferative SRF/Mkll target genes including mig6/ 
errfi-1, a negative regulator of the EGFR-MAPK pathway, 
were identified [65], or the tumor suppressor gene Eplin-a 
was described as a direct target of the SRF/Mkll path- 
way [66]. Furthermore, expression of a constitutively ac- 
tive form of Mkll in oncogenic ras- or src-transformed 
rat intestinal epithelial cells injected into the spleen of 
nude mice significantly suppressed tumor formation 
and reduced liver metastases by rescuing the expression 
of the SRF/Mkll targets tropomyosin and caldesmon 
[67]. In line with these findings, we could show that 
high expression of SRF/Mkll target genes is associated 
with an improved clinical outcome in breast cancer pa- 
tients. However, the opposite is the case for high expression 
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Figure 6 The SRF-independent/SAP-dependent genes represent a bad prognostic signature for breast cancer patients. Tumors in the 
1881 -sample breast cancer data set were stratified into tliree quantiles, low, intermediate and high, based on SRF-dependent/SAP-independent 
(A) or SRF-independent/SAP-dependent (B) gene expression. (A) Kaplan-Meier survival analysis using distant metastasis free survival (DMFS) as endpoint 
and 10-year censoring for all tumors (n = 1379; left panels), or in the subgroups of lymph node (LN)-negative (n = 1 1 1 1; middle panels) and untreated 
tumors (n = 821; right panels) was performed using the GOBO bioinformatics tool, interrogating the group of SRF-dependent/SAP-independent target 
genes. P-value is calculated using log-rank test. (B) Kaplan-Meier survival analysis for tumors with expression of SRF-independent/SAP-dependent Mkll 
target genes was performed as in (A). Association with clinical outcome was assessed in the subgroups of ER-positive (n = 856; left panel). 
Grade 1 (n = 1 41 , middle panel) and Grade 2 (n = 446; right panel) tumors in addition to all tumors and the subgroups used in (A). P-value is 
calculated using log-rank test. 



of SAP-dependent Mkll target genes. These genes are asso- 
ciated with poor clinical outcome predominantly in less ag- 
gressive tumors such as LN-negative, ER-positive, Grade 1 
and 2 tumors, which makes them valuable predictors of 
breast cancer progression. A scheme that depicts our model 



for Mkll action in breast cancer is presented in Figure 8. In 
this model Mkll is transactivating SRF-target genes in less 
aggressive tumors, while in the course of cancer progres- 
sion and metastatic behavior Mkll is activating a new 
group of genes in a SAP-dependent manner either by direct 
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Figure 7 Elevated expression of SAP-dependent Mkll target genes is a poor prognosis factor in breast cancer independent of histological 
grade. Multivariate analysis supporting the Kaplan-Maier survival analysis (shown in Figure 6) for the SRF-dependent/SAP-independent (A) and 
SRF-independent/SAP-dependent (B) gene sets, was performed using the GOBO bioinformatics tool. The analysis was executed for all tumors 
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only, or when not all of the investigated cases have clinical follow-up or clinical information for the specified covariate. 



interaction with the promoters of these genes or by inter- 
action with additional DNA-binding factors. 

Interestingly, in parental HCll cells many of the genes 
that we found in the SAP-dependent gene set that foster 
cell proliferation and migration and may cause poor 
survival of breast cancer patients are also induced by 
mechanical strain. A recent study has demonstrated 
that inhibition of cell spreading due to a lack of matrix 
stiffness is overcome by externally applied stretch, sug- 
gesting that similar mechanotransduction mechanisms 
sense stiffness and stretch [69]. Tumor stroma is typically 
stiffer than normal stroma. In breast cancer, diseased tis- 
sue can be 10 times stiffer than normal breast [70,71]. It is 
known that abnormal ECM stiffness plays an important 
role in cancer progression [72,73], but the mechanisms by 
which stiffness influences cancer progression are still 
under investigation. If we assume that we have discovered 
a general reaction of mammary epithelial cells to mechan- 
ical strain, we envisage that epithelial cells in a stiff. 



mechanically dynamic tumor environment may react by 
inducing a SAP-dependent Mkll gene set that in turn 
affects tumor progression. Furthermore, the products of 
these genes, many of which are involved in ECM turn- 
over and function, for example Lox [58], Mmps [74], 
Adamtsl6 [31] or Wispl [44] might themselves manipu- 
late the tumor microenvironment, thereby influencing 
tumor cell survival by a positive tumorigenic feedback 
loop. 

Finding how to switch the mode of action of Mkll be- 
tween SRF transactivation versus its SAP-dependent 
transcriptional activity is a subject of ongoing research 
in our lab that in future may help with the development 
of new therapeutic interventions for breast cancer. Post- 
translational modifications such as sumoylation are known 
to influence Mkll transcriptional activity [75] and phos- 
phorylation has been shown to influence interaction of 
Mkll with nuclear actin resulting in transcriptional 
changes [76,77]. Further characterization of these and 
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Figure 8 Schematic representation of the Midi action in breast cancer. A circular Mkll model is depicted with four of its domains: RPEL, 
actin binding motifs with RPxxxEL core consensus; Bl, basic domain involved in SRF-binding; SAP, homology domain found in the nuclear 
proteins SAF-A/B, Acinus, PIAS; TAD, transactivation domain. Serum response factor is drawn as a red shape and putative unidentified 
DNA-binding proteins as white shape with a question mark. Mkll exerts two distinct modes of action: one of them is through the Bl domain 
required for serum response factor (SRF)-binding activity and induction of SRF/Mkll target gene expression; the other one is strongly dependent 
on the SAP domain and triggers the expression of a specific set of pro-proliferative and pro-migratory genes that we called SAP-dependent Mkll 
target genes. High expression of SRF-dependent genes is associated with good clinical outcome for breast cancer patients, whereas elevated 
expression of SAP-dependent targets correlates with poor prognosis and indicates a significant role for these genes in breast cancer progression. 



other post-transcriptional changes of Mkll deserve spe- 
cial attention when trying to answer the above question. 

Conclusions 

In the current study, we discovered a breast cancer- 
specific set of genes that is highly interesting as a prog- 
nostic marker and therapeutic target for several reasons. 
(1) The expression of this gene set is regulated by Mkll 
and its SAP domain and is independent of SRF. (2) The 
SAP-dependent, SRF-independent Mkllsignaling is trig- 
gered by mechanical strain and may thus be activated in 
stiff tumors with a high stromal content and high inter- 
stitial tissue pressure. (3) This gene set is composed of 
interesting members some of which represent novel 
candidates for playing a functional role in cancer and 
others that have already been implicated in cancer- 
related functions, as for example tenascin-C, a meta- 
static niche component important for lung colonization 
[8], or Lox as a gene mediating collagen crosslinking re- 
sponsible for fibrosis-enhanced metastasis [58]. (4) The 
SAP-dependent Mkll target genes are associated with a 
poor clinical outcome in breast cancer patients, not re- 
ceiving adjuvant therapy or having a cancer classified 
as non-aggressive such as LN-negative, ER-positive, 
Grade 1 or 2 tumors. This makes these genes potential 
valuable prognostic markers in selecting patients who 



may benefit from an immediate and/or more aggressive 
therapy. 

Methods 

Cell culture 

Full length Mkll (FL-Mkll) and the two Mkll mutants, 
mutBl-Mkll comprising alanine substitutions of four 
amino acids (underlined) in the Bl domain of Mkll (KKA 
KELKPKVKKLKYHQYIPPDQKQD) [78] and ASAP-Mkll 
with a deletion of the SAP domain [15], were constructed 
based on transcript variant 1 (GenBank accession num- 
ber NM_153049) as previously described [13]. All Mkll 
variants were expressed as C-terminal RFP-tagged fusions. 
An empty vector expressing RFP alone was previously 
described [13]. 

HCll mammary epithelial cells, kindly provided by 
Dr. N. Hynes (Basel, Switzerland), were grown in RPMI- 
1640 medium supplemented with 10% FCS, 5 (ig/ml insu- 
lin (Sigma, Buchs, Switzerland) and 10 ng/ml epidermal 
growth factor (EGF; Invitrogen, Zug, Switzerland). In most 
of the experiments, the HCll cells were starved in 0.03% 
FCS/RPMI without EGF. To obtain HCll cells stably 
expressing FL-Mkll-RFP (HCll-FL), mutBl-Mkll-RFP 
(HCll-mutBl), ASAP-Mkll-RFP (HCll-ASAP) or RFP 
alone (HCll-empty vector), cells were transfected using 
FuGENE® 6 (Roche, Basel, Switzerland) and selected 
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with Geneticin (1 mg/ml; Roche) for 14 days before 
fluorescence-activated cell sorting (FACS) of RFP-positive 
cells on a Vantage SE (Becton Dickinson, Basel, Switzerland). 
Cell viability of the four HCll cell strains was assessed by 
the CellTiter-Blue viability assay (Promega, Duebendorf, 
Switzerland). 

Cell proliferation assay 

Proliferation rates of the HCll cell strains were determined 
using BrdU incorporation assay (Roche). After 24 h of star- 
vation, cells were plated in triplicate on Black 96-well mi- 
crotiter plates (PerkinElmer, Schwerzenbach, Switzerland) 
at 5 X 10^ cells/well in 3% FCS/RPMI and allowed to pro- 
liferate for 0, 24, 48, 72 and 96 h before labeling with BrdU 
for 2 h. BrdU incorporation into newly synthesized DNA 
was determined according to the manufacturer s protocol 
using a Luminometer Mithras LB940 (Berthold Technolo- 
gies, Regensdorf, Switzerland). Experimental values were 
normalized to the values of HCll- ASAP cells at the time 
point 0. Data represent means ± SD from three independ- 
ent experiments. 

Cell migration assay 

Cell migration was assayed using transwell polycarbonate 
membrane inserts (6.5 mm; Corning, Amsterdam, The 
Netherlands) with 8 [im pores as described [79]. After 
24 h of starvation, 5 x 10^ cells were plated in the top in- 
sert chamber with 100 \i\ serum-free RPMI. The lower 
chamber was filled with 600 \i\ 10% FCS/RPMI. Cells were 
allowed to migrate across the filter for 22 h at 37°C before 
fixation and crystal violet- staining. Images of duplicate in- 
serts were acquired on a Nikon Eclipse E600 using lOx 
magnification and a color CCD camera. Migration was 
quantified by measuring the area covered by migrated cells 
using the Fiji distribution of ImageJ [80]. Data represent 
means ± SD from three independent experiments. 

Mechanical stimulation of cells 

2 X 10^ HCll cells/well were seeded in BioFlex® 6-well 
culture plates (Flexcell International, Hillsborough, NC, 
USA) coated with either growth factor reduced-Matrigel 
(BD Biosciences, Basel, Switzerland) or fibronectin [11]. 
Cultures were starved for 24 h before applying either 
equibiaxial cyclic strain (15%, 0.3 Hz) or static strain 
(20%) at 37°C for 1 h using Flexcell FX-4000 (Flexcell 
International). Cells cultured under the same conditions 
and not exposed to strain were used as a resting control. 
After mechanical stimulation, cells were lysed and total 
RNA was isolated using the RNeasy Mini Kit (Qiagen, 
Basel, Switzerland). 

Transcript profiling and bioinformatics analysis 

HCll cell strains stably expressing Midi variants were 
starved for 48 h before total RNA was extracted. 



converted into labeled cDNA and hybridized to Affy- 
metrix GeneChip Mouse Gene 1.0 ST arrays. RMA- 
normalized expression values were calculated with the 
Affy package from Bioconductor 2.4 [81], and differen- 
tially expressed genes were identified using moderated 
^-statistics calculated with the empirical Bayes method 
as implemented in the Bioconductor limma package 
[82]. To be considered as differentially expressed be- 
tween HCll-FL and HCll-mutBl or HCll-ASAP 
cells, genes had to pass the filters: adjusted P-value < 
0.01 (with Benjamin-Hochberg false discovery correc- 
tion), a minimum absolute Unear fold change differ- 
ence of 2.0 and a minimum average expression value of 
4.0 (log2). Microarray data files are available from the 
Gene Expression Omnibus (GEO; http://www.ncbi.nlm. 
nih.gov/geo/), accession number GSE44907. Using the 
above parameters, gene lists of the two contrasts (mutBl/ 
FL and ASAP/FL) were compared resulting in the forma- 
tion of three gene groups: SRF-dependent/ SAP-independ- 
ent, SRF-dependent/SAP-dependent and SRF-independent/ 
SAP-dependent. The three gene sets were analyzed using 
the bioinformatics softwares: 1) IPA (Ingenuity® Systems; 
www.ingenuity.com); and 2) GOBO (http://co.bmc.lu.se/ 
gobo) [59]. In order to use the latter tool, Affymetrix Gene- 
Chip Mouse Gene 1.0 ST IDs were mapped to Aftymetrix 
Human Genome U133A IDs using Biomart for Ensembl 
build 66. The module "Gene Set Analysis Tumors" was 
used to investigate the expression pattern and to per- 
form survival and functional correlation analyses for the 
SRF-dependent/SAP-independent and SRF-independent/ 
SAP-dependent gene sets across 1881 breast cancers char- 
acterized by Aftymetrix Human Genome U133A arrays. 

RNA analyses by qRT-PCR 

Total RNA was isolated from HCll cell strains after 
24 h of incubation either in 0.03 or 3% FCS/RPMI. 
RNA was reverse transcribed and relative tenascin-C 
and c-fos mRNA levels were detected as described 
[12,13]. Relative mRNA levels for the genes listed in 
Table 1, normaUzed to Gapdh, were measured using 
Platinum" SYBR" Green qPCR SuperMix-UDG with 
ROX (Invitrogen) and the primers listed in Additional 
file 4: Table S4. Real-time PCR was performed in a Ste- 
pOnePlus Real-Time PCR System (Applied Biosystems, 
Rotkreuz, Switzerland) using a standard cycling profile. 
All samples were run in duplicate. Data were analyzed 
by the ACt method [83] and presented as fold changes 
in mRNA expression levels between HCll-FL and 
HCll-ASAP cells. RNA from stretched cells was ana- 
lyzed by qRT-PCR using the efficiency AACt method 
[84] that included a further normalization to the rest- 
ing control. Data represent means ± SD from three in- 
dependent experiments. 
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Additional file 4: Table S4. Primer sequences. Table S5. Promoter 
constructs. Figure SI. Quantification of SAP-dependent Mkll target gene 
expression using qRT-PCR analysis. Figure S2. Differential expression of 
newly discovered Mkll target genes in HCll strains overexpressing either 
FL-, mutBI- or ASAP-Mkll constructs (protein analyses performed by 
immunoblotting and zymography). Figure S3. SAP-dependent Mkll 
target genes are correlated with a very high significance (P < 0.00001) 
with the two proliferation modules - mitotic checkpoint and mitotic 
progression (a functional correlation analysis performed using the GOBO 
bioinformatics tool). 



Protein analyses by immunoblotting and zymography 

After 24 h of starvation, whole-cell extracts from the 
three HCll strains were prepared in RIPA buffer and 
immunoblotting was performed as described [12,13]. 
The following primary antibodies were used: mAb65F13 
anti-Mkll [12], MTnl2 anti-Tnc [85], anti-Wispl/CCN4 
(clone 214203, R&D Systems), anti-Nox4 (NBllO-58851, 
Novus Biologicals), anti-Vcl (clone hVIN-1, Sigma) and 
anti-Gapdh (ab9485. Abeam). 

After reaching 90% confluency, HCll strains were 
starved for 48 h before conditioned medium was col- 
lected, concentrated and analyzed by zymography as 
described [86]. 

Promoter-reporter assays 

The tenascin-C promoter used in this study was described 
as TNC 247 bp [13]. Promoters of Acta2 [87] and all SRF- 
independent/SAP-dependent genes described in Table 1 
were PCR-amplified from genomic DNA and corresponded 
to the sequences listed in Additional file 4: Table S5. Each 
promoter contained > 500 bp 5 ' of the TSS and was cloned 
into the pSEAP2-Basic (Clontech, Saint-Germain-en-Laye, 
France). For some promoters also 200 bp proximal pro- 
moter sequences were cloned as described above. All clones 
were verified by DNA sequencing. 

HCll cells in 6-well plates were cotransfected with 
1 [ig of the SEAP reporter vectors, 1 [ig of pcDNA3 vectors 
encoding Mkll variants [13], and 200 ng of the secreted 
luciferase MetLuc vector (Clontech) used to normalize 
for transfection efficiency. Cells were cultured in 0.03% 
FCS/RPMI for 24 h before enzymatic activity measure- 
ments were performed as described [13]. Experimental 
values represent averages of three independent experi- 
ments, each performed in duplicate. 

Statistical analysis 

Numerical results were expressed as means ± SD. Stat- 
istical analysis was completed using GraphPad InStat 
Software, version 3.05. The two-tailed Students t test 
was used to evaluate differences between two groups. 
Multiple comparisons were performed using one-way 
analysis of variance (ANOVA). Values of P less than 0.05 
were considered statistically significant. Statistics for bio- 
informatics analyses is given in figure legends. 

Additional files 



Additional file 1: Table 51. SRF-dependent/SAP-independent 
probeset list. 

Additional file 2: Table 52. SRF-dependent/SAP-dependent probeset 
list. 

Additional file 3: Table 53. SRF-independent/SAP-dependent 
probeset list. 
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